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We present the initial conditions for a decaying cosmological perturbation and study its signatures 
in the CMB anisotropies and matter power spectra. An adiabatic decaying mode in presence of com- 
ponents which are not described as perfect fluids (such as collisionless matter) decays slower than in 
a perfect-fluid dominated universe and displays super-Hubble oscillations. By including a correlated 
decaying mode with a red or a scale invariant spectrum, the anisotropy pattern shows super-imposed 
oscillations before the first Dopplear peak, while with a blue spectrum the amplitude of the sec- 
ondary peaks relative to the first one and the matter power spectrum can be altered. WMAP first 
year data constrain the decaying to growing ratio of scale invariant adiabatic fluctuations at the 
matter-radiation equality to less than 10%. 

PACS numbers: 98.80.-k, 98.80. Cq, 98.80.Es 



Introduction. The temperature and polarization 
anisotropy pattern of the Cosmic Microwave Background 
Radiation (CMB) is an invaluable snapshot of the initial 
conditions for the primordial perturbations responsible 
for structure formation. The most general initial con- 
dition in the linear regime can be written as a sum of 
several correlated modes that can be classified as adi- 
abatic or isocurvature depending on whether the total 
non- adiabatic pressure perturbation vanishes or not, and 
as growing or decaying depending on their time behavior. 
The standard paradigm states that the fluctuations be- 
gin in their pure adiabatic mode and that the decaying 
component becomes rapidly negligible. Any deviation 
from this simple one-mode description may be a signal 
of new phenomena in the early universe. For instance, 
additional isocurvature (growing) modes 0] may be the 
remnants of multi-field inflation [3, 0, 0| or of primordial 
phase transitions and a rich literature has developed on 
this topic. Although the most recent observations by the 
WMAP satellite have shown that the cosmological per- 
turbations start in an initial condition very close to the 
growing adiabatic one a sizable amount of isocurva- 
ture perturbations either in the neutrino or in the CDM 
component is still allowed by the data Q. 

In contrast, it is rather remarkable that, to the best 
of our knowledge, no attempt to constrain the decay- 
ing modes on the CMB or matter power spectrum has 
ever been performed; in fact, the CMB spectrum of the 
decaying mode has never been evaluated. All observa- 
tional constraints on the initial conditions neglected to 
take into consideration the adiabatic decaying mode, on 
the implicit assumption that they are completely negli- 
gible. There are however three main reasons why the 
decaying modes might be observable on the CMB. First, 
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if they are generated with a spectrum different from the 
growing ones, then they might be dominant either at 
small or large scales whatever their amplitude at a given 
scale is. Second, although the decaying modes die out 
during a cosmological epoch, any sharp transition to a 
different epoch characterized by different components or 
different equation of states will regenerate a new growing 
and a new decaying modes: the decaying modes trace 
then the occurrence of cosmological transitions. Third, 
if a non standard era occurred in the early universe (for 
instance a mix-master singularity or a bounce ja|), an adi- 
abatic decaying mode with a large initial amplitude may 
be present at the onset of the expanding stage. There- 
fore, it may be conceivable to have a non-vanishing decay- 
ing part when fluctuations reenter the Hubble radius in 
bouncing models. Aim of this paper is to discuss from a 
general and phenomenological point of view the general 
adiabatic mode and to derive constraints from a likeli- 
hood analysis of CMB data. In this new era of precision 
cosmology, this is a novel test for the coherence of pri- 
mordial cosmological fluctuations. 

The decaying initial conditions. In order to concretely 
discuss this subject it is useful to write the perturbed 
scalar part of Robertson- Walker metric both in the lon- 
gitudinal and in the synchronous gauge, respectively: 

ds 2 = a 2 (T)[-dT 2 {l + 2^) + dx l dx j ~i i] {l-2<$>)\ 

= a 2 (T)[-dT 2 + dx l dxH ll3 + hij)] , (1) 

where 7y is the flat background metric and a(r) is the 
scale factor (we fix the present value do = 1 and we 
have omitted the differences in the coordinates in the 
two gauges in Eq. (1)). The scalar metric perturbation 
in the synchronous gauge /ly is separated in its trace h 
and traceless potential 677. In a single perfect fluid dom- 
inated universe $ = 'J, and the Fourier components of 
the gravitational potential <& obey the following (driven) 
dissipative wave equation which admits two linearly in- 



2 



dependent solutions [7j: 
$" + 3H(l + w)$'-wV 2 $ 



[2H' + (1 + 3w)H 2 ] $ 



= 
(2) 

where w = p/p is the equation of state of the fluid and 
H is the conformal Hubble parameter. By considering as 
an example a radiation dominated universe (a(r) oc r), 
<I> admits as solution: 
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where y = kr/^/3. This means that by studying the 
impact of the decaying mode we are also testing the co- 
herence hypothesis of the primordial cosmological fluctu- 
ations, which is at the core of the inflationary scenario 
combined with a standard cosmology after nucleosynthe- 
sis. 

Let us now derive the initial conditions at large 
scales near or after the epoch of nucleosynthesis, around 
z ~ 10 s — 10 9 . At this stage baryons, CDM, photons 
are treated as perfect fluids, while neutrinos freestream 
since they already decoupled from the primordial plasma 
(around z ~ 10 11 ). We work from now on in synchronous 
gauge, as in most Boltzmann codes. The indexes b, c, r, v 
will denote baryons, CDM, photons and neutrinos, re- 
spectively. Without neutrinos the initial conditions to 
the lowest non trivial order are, up to an arbitrary over- 
all factor C [13]: 
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where 9 r is the peculiar velocity divergence. These equa- 
tions are to be complemented by S r — — 2h/3, S c = 5b = 
—h/2, 9 C = and 9b = 6 r , where D is the amplitude ratio 
decaying-to-growing and x = fcr (Note that during radia- 

— 1/2 

tion t w a e q' a/Ho if aeq is the equivalence scale factor, 
normalized in such way that do = 1 today. In principle 
this relation receives a small correction from the recent 
A-dominated epoch, but for semplicity we shall neglect 
it thereafter.). 

Adding neutrinos is not an easy task, since neutrino 
traceless pressure perturbations couple to the shear of the 
metric ($ ^ VP in the longitudinal gauge). The neutrino 
energy-density fraction R v is defined as 



ftv + P 1 



with p 7 , p„, are the energy density of neutrinos and pho- 
tons (after electron-positron annihilation) and N v is the 
number of neutrino species, respectively. On considering 
massless neutrinos R v is constant in time (this is also a 
very good approximation for realistic neutrinos masses) . 



The physical initial conditions in presence of neutrinos 
are (again up to an overall factor): 
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To these one must add S r — —2h/3, S c = Sb = —h/2, 
9 C = and 9b — 9 r as above. A similar result was also 
found by Zakharov [Sj for a system of radiation plus a 
collisionless component (see also 0): our solution include 
baryons and CDM as well. Notice that h, r\ and 9 r in © 
reduce to (@J when R v — > 0, up to a trivial rescaling. 
In the following we will refer to the growing modes as 
g-modes and to the decaying modes as c?-modes. We 
assume that C, D are multiplied by their initial power- 
law spectrum with powers n a , rid, respectively, while we 
take tp independent of k. This means that the two linearly 
independent decaying solutions have the same spectral 
index. 

Several comments are in order. First, the decaying ini- 
tial condition 10 is as adiabatic as one can get: only 5 V 
does not agree with the adiabatic relation 8^ — —2h/3 
since it is the only component which has non-vanishing 
anisotropic pressure terms. Second, the decaying mode 
contains two arbitrary constants, D and ip, instead of just 
one, as for isocurvature modes: this happens because the 
order of the system of differential equations is increased 
by using a first-order time equation (introduced by the 
Boltzmann hierarchy for the moments of the neutrino 
distribution function) for a v and the number of decay- 
ing modes is doubled. By considering for simplicity the 
system of radiation plus neutrinos, the order of the lin- 
ear differential equation for h (see Eq. (93) of a t 
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leading order in x is therefore increased from four to five, 
leading to the growing and two decaying modes presented 
here, plus the two gauge modes proper of the synchronous 
choice as can be seen from the master equation for h: 
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Third, the mode decays more slowly than in a perfect- 
fluid dominated universe (r -1 / 2 instead of t _1 in Eq. <@J 
in rf) and oscillates in log a; with frequency 7/2, when 
R v > 5/32. Oscillations are erased when R v < 5/32; 
however this does not occur here since we keep the conser- 
vative value N„ = 3. Fourth, the above initial condition 
imply that F u3 ~ O ((fca) 1 / 2 x cos, sin). Therefore the 
third and the higher moment of the neutrino distribution 
are regular in the past [F vm ~ O ((/ca) m ~ 5 / 2 x cos , sin) 
with m > 3) and can be safely neglected in the above 
initial conditions, as for the growing mode fltil |. 

Impact on observations. There are several novel effects 
due to the adiabatic decaying mode on the CMB temper- 
ature, polarization and matter power spectra. The oscil- 
lations in the gravitational potential will induce an oscil- 
lating pattern of SW and ISW effects at low multipoles 
(before the first peak). There is a chance that the waves 
at I < 200 seen in the first year WMAP temperature 
data are a signature of a non-vanishing decaying mode 
on large scales: we are exploring quantitatively this issue 
[Tlj . Note that the superimposed oscillations induced by 
a decaying mode appear only on scales larger the first 
acoustic peak, since the oscillations in log 2 found here 
are proper only of large scales. The prediction of super- 
imposed oscillations only before the first peak sharply 
distinguishes the effect of a decaying mode from design- 
ing oscillations on the primordial Fourier matter power 
spectrum |l2T |. 

The oscillations are present when g and d modes are 
correlated: the pure decaying spectrum, being quadratic 
in the oscillating functions, shows in fact attenuated d- 
oscillations. This effect will be maximized for rid < n g , 
because in this case the largest scales are enhanced with 
respect to the smaller ones, and also for small optical 
thickness r, because an early reionization erases part of 
the oscillation pattern. As is also clear from Fig. 1, the 
contribution from a correlated decaying adiabatic mode 
is in general not vanishing for £ larger than the first peak. 

Likelihood analysis. Multiplying the decaying mode by 
the decaying spectrum Pd ~ k nd , the fractional ampli- 
tude of mode k in the metric perturbation h can be writ- 
ten as Dfc(r) = L>(fc/A:*)" d / 2 (fcr)- 1 / 2 where we put k* = 
0.05/i/Mpc. We can perform therefore a quick estimation 
of the expected effect. For a flat spectrum (rid = 1) the 
fractional amplitude of the decaying mode with respect 
to the growing one can be approximated as D k * (r) = 
D(k*T)~^ 2 « DiHo^/a^/ak*) 1 / 2 » O-OlDa" 1 / 2 . Ex- 
trapolating to the equivalence time (aeq ~ 3x 10~ 4 (l3l |l 
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Figure 1: Temperature power spectra for D = 0.2 , ip = 
tt/4, n d = 1 (left) and D = 0.2 , ip = ,n d = 3 (right). 
The other parameters are n g = 1 and h = 0.7, fl^h 2 = 
0.025, Q c h 2 — 0.125. The solid line is the maximally cor- 
related spectrum, the dotted has no decaying mode, the 3 
dotted-dashed is the uncorrelated spectrum, the dashed is 
the decaying spectrum (multiplied by 3) and the dot-dashed 
the correlation. Note that in the left panel the correlation is 
negative on £ < 60 and positive on larger £. 
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Figure 2: Left and middle panel, TE and E power spectra for 
the same parameters of the left panel in Fig. 1 (the convention 
with lines are the same as in Fig. 1 and the pure decaying 
spectrum is not plotted). To the right, unnormalized power 
spectrum for CDM for a vanishing decaying mode (solid line), 
for na = 1 D = 0.2 , —0.2 (dotted and dashed, respectively), 
rid = 3, D = 0.2 , —0.2 (dot-dashed and 3-dot-dashed). 



, this amounts to Dk> (aeq) ~ 0.56 D. As from Fig. 1, 
values of D k * (a e q) ~ O(10 _1 - 10~ 2 ) (with rid ~ n g w 1) 
lead to waves super-imposed to the power spectrum of a 
few percent amplitude, similar to those seen in the first 
year WMAP data. This also shows that for rid > n g the 
decaying mode may lead to effects on small scales, as can 
be seen from Fig. 2, and can be constrained by Planck 
Q and LSS data. 

In a subsequent paper we will investigate in detail the 
full likelihood of the model including the features of the 
low multipole waves In this paper we present only 
the main result: the constraints from the WMAP data 
on the decaying amplitude D. To simplify the compari- 
son to real data we assume here that the g — d modes are 
either maximally correlated, uncorrelated or maximally 
anti-correlated . There is no loss of generality by as- 
suming therefore D > 0. Inserting the initial conditions 
(jHJ in the Boltzmann code produces maximally corre- 
lated spectra; the uncorrelated and anticorrelated cases 
are obtained by linear combinations of a basis of three 
maximally correlated srjectra (for each combination of 
parameters), as in e.g. 
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Figure 3: Likelihood for the decaying amplitude D for various 
choices of rid and of the g — d correlation as in the labels, 
marginalized over the other parameters. 

Values of rid substantially larger than unity generate 
strong deviations at small scales, both in the Ci and in 
the matter power spectrum, so they are expected to be 
severely constrained by large scale structure data and by 
large-^ data [ijj. Fig. 2 shows for instance that the mat- 
ter power spectrum for rid > 3 is extremely deficient at 
small scales, so we exclude such values from the analy- 
sis. We limit ourselves to discuss a few cases: a scale 
invariant rid — 1 and two blue spectra, rid = 1.5,2. We 
explored a grid of over 10 6 models spanning eight param- 
eters (for each choice of the d — g correlation): D,rid,ip 
plus the standard parameters n s , h,Q(,h 2 ,Q c h 2 , t (the 
overall amplitude is automatically integrated out). We 
restrict ourselves to the massless neutrinos case. We 
compare the predictions obtained from a modified ver- 
sion of CMBFAST 0| with the first year temperature 
and polarization WMAP data (Ta j. The global best 
fit is for D — I, rid — 2, tp = —0.785 rad , no g — d 
correlation, and n s = 0.98, h = 0.69, tt b h 2 = 0.024, 
n c h 2 = 0.127, t = 0.15, with X 2 = 970 (for the tempera- 
ture spectrum alone) which is however just marginally 
better (Ax 2 = 4) than the standard ACDM best fit. 
The results of the likelihood analysis are shown in Figs 
3. For a flat spectrum (rid — 1) the data constrain 



D < 0.2 at 95% CL both for correlated and anticorre- 
lated modes and D < 0.1 for uncorrelated modes. This 
means that any process that excited a decaying mode 
of the same amplitude and spectrum of the growing one 
(i.e. Dfe*(a«) = 1) has occurred not later than a*, where 
a* = 10~ 4 D 2 < 4 • 10~ 6 . For larger rid the upper limit to 
D moves to larger values; for rid = 2, an amplitude D of 
order unity is allowed in the uncorrelated case. As it can 
be seen from Fig. 3, the effect of the g — d correlation 
is rather limited. Similarly, within the range of values of 
rid we adopted here the phase ip has a very small effect 
on the likelihood (a larger effect is obtained for rid < !)• 
Finally, the estimation of the other cosmological param- 
eters is not much affected by the addition of d-modes 
[Tlj . Therefore, it is safe to consider D < 1 as an upper 
limit independent of all the other parameters involved in 
the model. Fixing rid = 1 the upper limit strenghtens to 
L><0.2. 

Conclusions. We have derived the initial conditions 
for an adiabatic decaying mode in a standard primordial 
soup (neutrinos, baryons, CDM and radiation) in a spa- 
tially flat universe in Eq. . We have shown that decay- 
ing modes induce super-Hubble oscillations and observ- 
able deviations at small scales in the CMB and matter 
power spectra, depending on their amplitude and spec- 
trum. WMAP first year data alone constrain D < 0.2 at 
95% CL for a scale invariant spectrum almost regardless 
of the other parameters. This means that scale invari- 
ant adiabatic fluctuations should consist of their grow- 
ing part for at least ~ 90% at matter-radiation equality. 
When rid > n g , larger values for D are still allowed by 
observations. 
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